Instabilities of wormholes and regular black holes supported by a phantom scalar field 
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We test the stability of various wormholes and black holes supported by a scalar field with a 
negative kinetic term. The general axial perturbations and the monopole type of polar perturbations 
are considered in the linear approximation. Two classes of objects are considered: (i) wormholes 
with flat asymptotic behavior at one end and AdS on the other (M-AdS wormholes) and (ii) regular 
black holes with asymptotically de Sitter expansion far beyond the horizon (the so-called black 
universes). A difficulty in such stability studies is that the effective potential for perturbations 
forms an infinite wall at throats, if any. Its regularization is in general possible only by numerical 
methods, and such a method is suggested in a general form and used in the present paper. As a 
result, we have shown that all configurations under study are unstable under spherically symmetric 
perturbations, except for a special class of black universes where the event horizon coincides with 
the minimum of the area function. For this stable family, the frequencies of quasinormal modes of 
axial perturbations are calculated. 
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I. INTRODUCTION 

Modern observations [1] indicate that the Universe is 
expanding with acceleration. The most favored explana- 
tion of this acceleration is nowadays that the Universe 
is dominated (to about 70 %) by some unknown form of 
energy density with large negative pressure, termed dark 
energy (DE), while the remaining 30 % consisting of bary- 
onic and nonbaryonic visible and dark matter. It is often 
admitted that DE can be modeled by a self-interacting 
scalar field with a potential. Such a field acts as a neg- 
ative pressure source; it is called quintessence if its pres- 
sure to density ratio p / p — w > —1 and a phantom field 
if w < —1 while w ~ —1 corresponds to a cosmologi- 
cal constant A. Since observations admit a range of w 
including w = —1, all sorts of models are under study. 

One should note that values w < — 1 seem to be not 
only admissible but even preferable for describing an in- 
creasing acceleration, as follows from the most recent es- 
timates: w — —1.10 ± 0.14 (Icr) [2] (according to the 
7-year WMAP data) and w = -l-OQ^tl ^^l [3] (mainly 
from data on type la supernovae from the SNLS3 sam- 
ple) . In this connection, cosmological models with phan- 
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tom scalar fields, i.e., those with a negative kinetic term, 
have gained considerable attention in the recent years [4] . 

If such fields can play an important role in cosmology, 
it is natural to expect that they manifest themselves in 
local phenomena, for instance, in the existence and prop- 
erties of black holes and wormholes [5] . Quite a number 
of scalar field configurations of this kind have been de- 
scribed in the literature, see, e.g., examples of black holes 
with scalar fields (the so-called scalar hair) in [6] and 
wormholes supported by scalar fields in [7-9] and refer- 
ences therein. Thus, in [9] it was shown that, in addition 
to wormholes, a phantom scalar can support a regular 
black hole where a possible explorer, after crossing the 
event horizon, gets into an expanding universe instead of 
a singularity. Thus such hypothetic configurations com- 
bine the properties of a wormhole (absence of a center, a 
regular minimum of the area function) and a black hole 
(a Killing horizon separating R and T regions). More- 
over, the Kantowski-Sachs cosmology that occurs in the 
T region is asymptotically isotropic and approaches a de 
Sitter regime of expansion, which makes such models po- 
tentially viable as models of our accelerating Universe. 
Such configurations, termed black universes, were later 
shown to exist with other scalar field sources exhibit- 
ing a phantom behavior such as k-essence [10] and some 
brane world models [11] (in the latter case, even without 
a phantom field) . 

Both wormholes and black universes have been shown 
to exist as well in models where a scalar field exhibits 
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phantom properties only in a strong-field region while in 
the weak-field region it obeys the canonical field equation 
(the so-called trapped- ghost models) [12, 13] 

To sec whether or not such solutions can lead to vi- 
able models of black holes and wormholes, one needs to 
check their stability against various perturbations. Pre- 
viously, gravitational stability as well as passage of radi- 
ation through wormholes supported by a phantom scalar 
field were considered, in particular, in [14-18], with a spe- 
cial emphasis on massless wormholes (see also references 
therein) . 

In the present paper, we consider the linear stability 
of various static, spherically symmetric solutions to the 
field equations of general relativity with minimally cou- 
pled scalar fields, describing compact objects of interest 
such as asymptotically flat and AdS wormholes (M-AdS 
wormholes, for short, where M stands for "Minkowski") 
and black universes (in other words, M-dS regular black 
holes), and use as examples solutions obtained in [9]. 
We show that M-AdS wormholes are unstable in the 
whole range of the parameters while among black uni- 
verses there is a stable subfamily which corresponds to 
the event horizon located precisely at the minimum of 
the area function. 

The particular solutions whose stability is studied here 
arc certainly not general. A more comprehensive study 
is prevented by the fact that a sufSciently general solu- 
tion describing self-gravitating scalar fields with nonzero 
potentials is unknown, therefore it seems to be a nat- 
ural decision to study the properties of known special 
solutions. On the other hand, in [9] all possible reg- 
ular static, spherically symmetric solutions to the field 
equations were classified for phantom minimally coupled 
scalar fields with arbitrary potentials. One can see that 
the solution studied here, being certainly special, still 
represents a very simple but quite typical example repro- 
ducing the generic features of such configurations with 
flat, dS and AdS asymptotics; its another advantage is 
that it reproduces as a cpecial case the well-known El- 
lis wormhole, for which the stability results are already 
known [14-16]. 

The paper is organized as follows. Sec. II presents the 
backgrounds to be considered. Sec. Ill develops a general 
formalism for axial gravitational and Maxwell field per- 
turbations in a static, spherically symmetric background. 
Sec. IV is devoted to polar spherically symmetric pertur- 
bations. Sec. V discusses the stability of the black uni- 
verses and wormholes under consideration and analyzes 
the quasinormal radiation spectrum for the cases where 
the background configuration is linearly stable. In addi- 
tion, we there develop a numerical tool for reducing the 
wavelike equation with a singular potential to the one 
with a regular potential. In Sec. VI we summarize the 
results and mention some open problems. 



II. STATIC BACKGROUND 
CONFIGURATIONS 

Let us consider Lagrangians of the form 

L = ^g{R + e5«'5</.;„</.;^ - 2F(</.) - F^^F'"'), (1) 

which includes a scalar field, in general, with some po- 
tential V{4)), and an electromagnetic field F^,^; e = ±1 
distinguishes normal, canonical scalar fields (e = +1) and 
phantom fields (e = —1). In what follows, we present 
a perturbation analysis for static, spherically symmetric 
solutions for this general type of Lagrangian and then 
study the stability of some particular (electrically neu- 
tral) solutions. 

The general static, spherically symmetric metric can 
be written in the form 



ds"^ = A{r)dt'^ - A{r)-^dr'^ - R{rfdn'^ 



(2) 



where dO^ = d6'^ + sin^ 9dip'^ is the linear element on a 
unit sphere. 

Wc will consider the following static background [9]: 
i?(r) = (r^ 6^)1/2, 6 = const >0, 



A{r) = (r^ + 62) 



c 1 

62 + 62 + ^2 



-I- 



3m 



br 



63 \62-Fr2 



+ tan 



= \/2tan-i(r/6). 



(3) 



It is a solution to the Einstein-scalar equations that fol- 
low from (1) with F^^, = and the potential 



3to 



c 
'62 



2cos^ 



V2 



3 sm —j= cos —j= H — -7= 13 — 2 cos —= 
\2 \/2 v2 \ v2 



• (4) 



The solution behavior is controlled by the scale 6 
and two integration constants: c that moves the curve 
Blr) = A/ up and down, and ra showing the position 
of the maximum of B{r). Both R{r) and B(r) are even 
functions if m = 0, otherwise B{r) loses this symmetry. 
Asymptotic flatness at r = -|-oo implies 



26c: 



(5) 



where m is the Schwarzschild mass defined in the usual 
way. 

Under this asymptotic flatness assumption, for m = 

c = 0, we obtain the simplest symmetric configura- 
tion, the Ellis wormhole [7]: A = 1, V = 0. With 
m < 0, we obtain a wormhole with an AdS metric at 
the far end, corresponding to the cosmological constant 
^(0)| _ = V- < 0. Further on, such configurations 
will be referred to as M-AdS wormholes (where M stands 
for Minkowski, the flat asymptotic). Assuming m > 0, 
at large negative r we obtain negative A{r), such that 
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\A{r)\ ^ R^{r), and a potential tending to V- > 0. Thus 
it is a regular black hole with a de Sitter asymptotic be- 
havior far beyond the horizon, precisely corresponding to 
the above description of a black universe. 

In black- universe solutions, the horizon radius depends 

on both parameters m and b (min i?(r) = 6), which also 
plays the role of a scalar charge since 0/-\/2 «7r/2 — 6/r 
at large r. Since ^(0) = 1 + c, the minimum of R{r), 
located at r = 0, occurs in the R region if c > —1, i.e., 
if Sttto < 2b (it is then a throat, like that in wormholes), 
right at the horizon if c = — 1 (i.e., 37rm = 2b) and in the 
T region beyond it if c < —1, that is, 37rm > 2b. It is then 
not a throat, since r is a time coordinate, but a bounce 
in one of the scale factors R{r) of the Kantowski-Sachs 
cosmology; the other scale factor is A{r). 

Let us mention that another important case of the sys- 
tem (1), the one with V = 0, has already been stud- 
ied in a number of papers. In this case we are dealing 
with Fisher's famous solution ([19], 1948) for a canon- 
ical massless scalar in general relativity (e = and 
three branches of its counterpart for e = — 1, sometimes 
called the anti-Fisher solution, found for the first time 
by Bergmann and Leipnik [20] in 1957 and repeatedly 
rediscovered afterwards (as well as Fisher's solution). In 
the latter case the solution consists of three branches, one 
representing wormholes [7, 8] , the other two also contain- 
ing throats but with singularities or horizons of infinite 
area at the far end instead of another spatial infinity 
[8, 21]. An instability of Fisher's solution under spheri- 
cally symmetric perturbations was established long ago 
[22], a similar instability of the wormhole branch was 
discovered in [14, 16] and the same for the other two 
branches in [23] . The case of zero mass in the anti-Fisher 
solution represents the Ellis wormhole and is common 
with the solution under study, (3), m = 0. 

In what follows, we will assume m and use the 
mass m (which has the dimension of length and is equal 
to half the Schwarzschild radius in the units employed) as 
a natural length scale, putting, for convenience, ]3m] = 1. 
Then the constant c is used as a family parameter while 
b is found from the relation (5). 

We should remark that in all particular examples to be 
tested here for stability the background electromagnetic 
field is zero, but the general perturbation formalism is 
developed in Sec III for axial perturbations of systems 
with non-zero i^^^ as well. Actually in Sec IV non-zero 
is also allowed, simply because the monopole per- 
turbations do not excite an electromagnetic field in a 
spherically symmetric background. Perturbations of the 
electromagnetic field appear in higher multipoles of po- 
lar modes. The general formalism developed for non-zero 
Ff^,y is intended to be used in further studies of systems 
with both scalar and electromagnetic fields. 



III. LINEAR AXIAL PERTURBATIONS: 
GENERAL ANALYSIS 

In our analysis of axial perturbations, we use Chan- 
drasekhar's notations: cc" = t, ~ 0, = r, x^ = 9, 
so that the coordinates along which the background has 
Killing vectors are enumerated first. Following Chan- 
drasekhar [24], we consider the metric (2) as a special 
case of the metric 

ds^ = e^'^dt'^ - e2'''(# - adt - q2dr - q^dOf 

- e'^i'^dr'^ - e^^'de^. (6) 

Thus in (2) 

e'^-'=Air), 6^"= = A-i(r), 

The background electromagnetic field is taken in the 
form 

Fo2 = -Q*/R{rf, (8) 

that is, only a radial electric field, and is the (effective) 
charge. 

All calculations and results can be easily rewritten for 
magnetic fields, with nonzero Fl3, owing to the Maxwell 
field duality. One can bear in mind that configurations 
like wormholes and black universes can possess electric or 
magnetic fields without any real electric charges or mag- 
netic monopoles, due to their geometry, actually realizing 
Wheeler's concept of a "charge without charge" . 

It is easy to see that axial perturbations of a scalar 
field vanish. Then, ct, q2 and gs are perturbed, while V, 
112, IJ'S and u remain unperturbed. 

The axial gravitational perturbations obey the equa- 
tions 

SRi3 = 0, SRi2 = 2Q,R-^Foi (9) 
Let us introduce new variables: 

Qik = qi,k - (lk,i, QiO = Qifi - O",*, (10) 

with i,k = 2, 3, and 

£ = Foisin6'. (11) 

Recall that we use the numbers 0, 1,2,3 for t, cj), r and 
9 coordinates, respectively. The Maxwell equations, sub- 
ject to only first-order perturbations, have the form 

{e^+i'^F^2),3 + {e^+'''F3r),2 = 0, (12) 

(e^+''Foi),2 + (e^+''^Fi2),o = 0, (13) 

(e^+'^Foi),3 + (e'^+^^Fi3),o = 0, (14) 

(e'^^+'^^J^oi),o + (e"+^-^^^i2),2 

+ (e'^+''^Fi3),3 = e^+''^i^o2Qo2. (15) 
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After some algebra the Maxwell equations can be writ- 
ten in the form 



Re-^E^Ofi - {e^''{Re''E),r),r + — sin^ E^e- 



1 



R \ ' sine, 
= -Q* ((7,2,0 - 92,0,0) sin^ & (16) 
Prom SRi3 = and SR12 = 2Q^R-^Foi it follows 

and 

= e^^+''+f''Q,R-^Foi. (18) 
After introducing the new function 

Q = iJ^AQas sin^ 6*, (19) 
Eqs. (17) and (18) can be written in the form 
A dQ 



i?2 sin^ e dr 
A dQ 



O-,3,0 - 93,0,0, (20) 



"O',2,0 + 92,0,0 



i?2 sin^ 



(21) 



Let us differentiate Eq. (20) in r and Eq. (21) in 6 and 
then add the results. After some algebra we have 



= 4Q„e''i?sin'^6' 



d / E 



de V sin^ 



(22) 



Now let us return to the Maxwell field perturbation 
equation (16). Using (21), we can get rid of the term 
containing 0-,2,o — 92,0,0 in (10)- After some algebra and 
using the relations E ~ e*"*, Q ~ e*'^*, we obtain 



{e^''iRe''E)^r),r + Eiuj^Re-" - AQle^R-^^ 



- e^R-^sine { E 



smt 



Q* dQ 
RMnO'de 



■ (23) 



The angular variable can be separated by the following 
ansatz: 



Q{r,0) = Q{r)C-^l\e), 



(24) 



= ^ M = mr)CiU{Ol (25) 



where are Gegenbauer polynomials. Eqs. (22) and 
(23) then read 

d/AdQx 2^n^^n foa\ 

{e^%Re''E)^r),r - (m^ + 2)e''R-^E 

+ Eiuj^Re-" - AQle^R-^) = -Q^QR-\ (27) 



where A = R^e^", /i^ = - 1)(£ + 2), and we use the 
tortoise coordinate r* defined by d/dr* = AR~^d/dr. 

After passing over from Q and E to the new functions 
Hi and H2 using the relations 



Q = RH2, Re^E = 
Eqs. (26) and (27) can be reduced to 



El 
'211' 



(28) 



R 



R^ 



where we have introduced the operator 



(29) 

(30) 

(31) 



For an electrically neutral background Q* — 0, Hi = 0, 
and Eq. (29) reduces to the Schrodinger-like equation 



fH2 

drl 



+ (w2 _ v^s{r))H2 = 



with the effective potential 



Kff(r) = A{r)- 



i?2 



+ R{R ^),r.,r, ; 



(32) 



(33) 



As a partial verification of the above relations, we can 
check that Eqs. (29) and (30) reproduce some known spe- 
cial cases. Thus, Eq. (33), obtained here in the Chan- 
drasekhar approach [24], coincides with Eq. (4.13) of [25], 
obtained in the Regge- Wheeler approach. In addition, 
Eqs. (29), (30) coincide with Eqs. (144), (145), p. 230 of 
[24] in the Rcissner-Nordstrdm limit, and, as a result, 
the potential (33) coincides with the well-known Regge- 
Wheeler potential for Schwarzschild black holes. 



IV. POLAR PERTURBATIONS 

For polar perturbations let us consider the Einstein- 
scalar equations following from (1) with F^i, = 

Rt^u + ea^$a,$ - 6f,,V{(^) = 0, 

en$ + V^=0, V^ = dV/dc/), (34) 

where □ = V"Vq is the d'Alembert operator. In the 

polar case we can put <t = (72 = 93 = 0, while Su, (5/i2, 
(5/U3, and Sip do not vanish. In addition, we perturb the 
scalar field, 
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Let us restrict ourselves to the lowest frequency modes, 
corresponding to spherically symmetric (or radial) per- 
turbations, £ = 0. This gives Stp = S^s. Then we can use 
the gauge freedom and fix^ 

6ip = 6fj,3 = =^ SR{r, t) = 0. (35) 

Tedious calculations allow us to reduce the perturba- 
tion equations to a single wave equation 

g-2/.2-2..j^ _ _ ^^/(^/ _ + + U5(j) = 0, (36) 
where 



Introducing the new function ^ by putting 

5(j, = ^fe^^+i"*, (38) 

we bring the wave equation to the Schrodinger-like form 



^ + (a;2 - V,s{n))^ = 0, 



with the effective potential 

rr 1 
VeS = U + 



R drl ■ 



(39) 



(40) 



^^(i^V-l).^.eW (41) 



V. STABILITY ANALYSIS 

A. Methods 

Finite potentials. If the effective potential 14ff is finite 
and positive-definite, the differential operator 



W 



(42) 



is a positive self-adjoint operator in Li^ {r , dr ) , the space 
of functions satisfying proper boundary conditions. Then 
W has no negative eigenvalues, in other words, there are 
no normalizable solutions to the corresponding Schrodin- 
ger equation with wcll-behavcd initial data (smooth data 
on a compact support) that grow with time, and the sys- 
tem under study is then stable under this particular form 



^ See more details on gauges and gauge-invariant perturbations in 
[14, 23]. The present notations are related to those in [23] as 
follows: 

r ^ u, fc" 1-^ 7, R{r) i-> r(u), (j,2 i-^ a, /is i-^ /3. 



of perturbations. Therefore, in our stability studies, our 
main concern will be regions where the effective poten- 
tials are negative since possible instabilities are indicated 
by such regions. 

The response of a stable black hole or wormhole to ex- 
ternal perturbations is dominated at late times by a set of 
damped oscillations, called quasinormal modes (QNMs). 
Quasinormal frequencies do not depend on the way of 
their excitation but are completely determined by the 
parameters of the configuration itself. Thus quasinormal 
modes form a characteristic spectrum of proper oscilla- 
tions of a black hole or a wormhole and could be called 
their "fingerprints". Apart from black hole physics, 
QNMs are studied in such areas as gauge/gravity cor- 
respondence, gravitational wave astronomy [26, 27] and 
cosmology [28]. 

The quasinormal boundary conditions for black hole 
perturbations imply pure incoming waves at the horizon 
and pure outgoing waves at spatial infinity. For asymp- 
totically flat solutions, the quasinormal boundary condi- 
tions are 



±oo. 



(43) 



The proper oscillation frequencies correspond in a 
sense to a "momentary perturbation", that is, to the sit- 
uation where one looks at a response of the system (say, 
a wormhole) to initial perturbation when the source of 
the perturbation stopped to act. This is the essence of 
the word "proper" . Thus, in the wormhole case, no in- 
coming waves are allowed coming from cither of the in- 
finities. The issue of the boundary conditions for quasi- 
normal modes of wormholes is not completely new and 
was considered in [29]. 

Therefore, for wormholes the condition "pure incom- 
ing waves at the horizon" is replaced with "pure outgoing 
waves at the other spatial infinity" [29, 30]. For asymp- 
totically anti-de Sitter black holes or wormholes, the AdS 
boundary creates an effective confining box [26], so that 
on the AdS boundary one usually requires the Dirichlet 
boundary conditions 



* 0, 



r* DO. 



(44) 



This choice is not only dictated by the asymptotic of 
the wave equation at infinity but is also consistent with 
the limit of purely AdS space-time: the QNMs of an 
AdS black hole approach the normal modes of empty 
AdS space-time in the limit of a vanishing black hole 
radius [31]. Thus quasinormal modes of a compact object 
(a black hole or wormhole) in AdS space-time look like 
normal modes of empty AdS space-time "perturbed" by 
the compact object. 

If the effective potential is negative in some region, 
growing quasinormal modes can appear in the spectrum, 
indicating an instability of the system under such pertur- 
bations. It turns out that some potentials with a negative 
gap still do not imply instability. If there are no growing 
quasinormal modes in the black-hole or wormhole spec- 
trum, this object is stable against linear perturbations. 
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Singular potentials and their regularization. It 

can be seen from (40) (provided that VR^ < 1) that 
the effective potential Ves for radial perturbations forms 
an infinite wall located at a throat, with the generic 
behavior Ves ~ (r — rthroat)"^, since we have there 
R'{r) ~ r — rthroat- As a result, perturbations are in- 
dependent at different sides of the throat, necessarily 
turn to zero at the throat itself, and we thus lose part 
of the information on their possible properties. To re- 
move the divergence, one can use a method on the basis 
of S-deformations as described in [32, 33] for solutions 
with arbitrary potentials V{(j)). For anti-Fisher worm- 
holes {V{(j)) = 0) it was applied in [14]. 

If the scalar potential V{(f)) is zero, the effective poten- 
tial (40) takes the form 



VMr) = ^ 



i?2 



eAf^ 



(45) 



In this case the general static solution (w = 0) to Eq. (39) 
is 



*o(r-) = CiR{l- 



R 



A'R' 
eARcj)'^ 
A'^^R' 



2A 



(46) 



and its special cases with Ci = and C2 = were used 

in [14, 23] to remove the divergence in 14fj. Indeed, we 
can introduce the new wave function 

* = 5* - — , where S = — -r-^, (47) 

which satisfies the equation 



+ - ^eff(r.)) * = 
with the effective potential 

Veff = 25^ — VeS, 



(48) 



(49) 



and the new effective potential is everywhere regular if 
either Ci = or C2 = 0. This transformation was used 
to prove the instability of anti-Fisher wormholes [14] and 
other anti-Fisher solutions [23], 

It turns out that C\ and C2 can be both nonzero, 
leading to regularized effective potentials with the same 
quasinormal spectrum, though not preserving the sym- 
metry (p -H- —(p. This means that we can fix the boundary 
conditions arbitrarily and, if we find a "good" static so- 
lution, we can use it to remove the singularity of the 
effective potential. 

For nonzero scalar potentials V{(p), analytical expres- 
sions for static perturbations are unknown. Therefore, 
such suitable solutions must be found numerically under 
proper boundary conditions. 

First we notice that since ^0 is a static solution to 
Eq. (39), S satisfies the Riccati equation 



dS_ 
dr* 



Substituting (40) into (50), we find an expansion for S 
near the throat as [23] 



S{r) 



1 + c 



4c2 4c2(4c2 + (l + c)7r2) 
H — TTTT-^ — ^ — — -r 



3m7r2 ' 9m2(l-hc)7r4 



(51) 



+ 52 - Feff = 0. 



(50) 



where K is an arbitrary constant. 

With this expression we find the boundary condition 
for the function S{r) at some points close to the throat on 
both sides. Then we integrate (50) numerically and find 
5(r) between the throat and both asymptotical regions. 
Having S{r) at hand, we find the regularized effective 
potential (49); it is finite at the throat due to (51). 

We assigned different values to the free constant K. 
As a rule, we were able to integrate Eq. (50) numerically 
in a sufficiently wide range of r near the throat. The 
regularized effective potentials obtained in this way al- 
ways lead to the same growth rate of the perturbation 
function. However, for some particular values of the 
numerical integration with the Wolfram Mathematica® 
built-in function encounters a growing numerical error. 
Apparently, in these cases some alternative methods of 
integration should be used. 

In addition, we used the same method to find numer- 
ically the regularized potential for the Branch B anti- 
Fisher solution whose stability properties had been stud- 
ied previously [23]. Although the numerically found po- 
tential differs from Eq. (68) of [23] , the time-domain pro- 
file again shows the same growth rate. This confirms the 
correctness of the method suggested and used here. 



B. M-AdS wormholes with negative mass 



Polar perturbations. Now we are in a position to 
apply the above method to the special case of M-AdS 
wormholes (3) with m < 0, which are asymptotically 
AdS as r — >■ — 00. In Fig. 1 one can see that although 
the initial (singular) cffcictivc potentials VeS are positive 
definite, the regularized potentials V^s are negative in 
some range around the throat and at the AdS boundary. 
This behavior usually indicates an instability but does 
not guarantee it [34]. Therefore, to prove the instability 
of M-AdS wormholes, we have used the time-domain in- 
tegration method proposed by Gundlach ct al. [35] and 
later used by a number of authors (sec, e.g., [36]). The 
method shows a convergence of the time-domain profile 
with diminishing the integration grid and increasing the 
precision of all computations. We imposed the Dirichlet 
boundary conditions at the AdS boundary, as described 
in [37]. Fig. 2 shows examples of time-domain profiles 
for the evolution of perturbations. The growth of the 
signal allows us to conclude that such wormholes arc un- 
stable against radial perturbations. At larger c we found 
the regularized potential with deeper negative wells and 
observed a quicker growth of the signal. Therefore we 
conclude that all such M-AdS wormholes are unstable. 
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FIG. 1: Effective potentials for spherically symmetric perturbations of M-AdS wormholes with m = —1/3, c — 0.3 (red), c — 0.5 
(green), c — 0.8 (blue) with divergences at the throat (left figure) and the corresponding regular ones found numerically (right 
figure). Larger values of c correspond to larger absolute values of Veff at the AdS boundary (they are positive for divergent 
potentials and negative for regular ones). 




FIG. 2: Time-domain evolution of spherically symmetric per- 
turbations of M-AdS wormholes with m — —1/3, c — 0.5 
(green, bottom), and c = 0.8 (red, top). 



Axial perturbations. The effective potentials Vcs for 
axial perturbations of M-AdS wormholes Vcs are plotted 
in Fig. 3. One can see that above some threshold value 
of c, Vcs{r) has a negative gap. This threshold value of 
c is c w 1.737 in units for which m = —1/3. A further 
increase of c makes the negative gap deeper, however, 
time-domain profiles for the evolution of axial perturba- 
tions show a decay of the signal without any indication 
of instability (see Fig. 4). 



C. Black universes 

Polar perturbations. Black universes correspond to 
the metric (2), (3) with m > 0, c < 0. Black universes 
with c < —1, or equivalently Snm > 26, have no throat 



in the static region. From Fig. 5 one can see that for 
c = — 1 the effective potential has a negative gap, how- 
ever, time-domain integration proves that in this case the 
black universe is stable. For smaller values of c, an ad- 
ditional negative gap appears between the peak and the 
horizon, leading to an instability even for c — —1.001. 
For large negative c the potential peak disappears, and 
the potential becomes negative everywhere (see Fig. 6), 
which inevitably creates an instability. 

At c > — 1 the throat r = is in the static region. In 
this case we numerically find the regularized effective po- 
tentials which have large negative gaps leading to insta- 
bilities (see Figs. 7, 8). As c approaches zero, the growth 
rate of time-domain profile decreases, still remaining pos- 
itive because the effective potential remains negative in 
a wide region near the throat (Fig. 7). 

As |c| grows the negative gap becomes deeper and nar- 
rower, giving way to a small positive hill, which becomes 
broader (see Fig. 8) as c approaches —1. However, we 
do not observe a decrease in the growth rate as expected 
when approaching the parametric region of stability near 
c = — 1. This can be an indication that the only param- 
eter for which a black universe can be stable is c = —1. 

Axial perturbations and quasinormal oscillation 
frequencies. From Fig. 9 we can see that in the static 
region the effective potential is positive-definite. There- 
fore, if we perturb a black hole "on the right" of the event 
horizon (Fig. 9), such perturbations are stable. Beyond 
the event horizon, in the cosmological region, the effec- 
tive potential can take large negative values, but this has 
no effect on perturbations propagating outside the hori- 
zon. We therefore conclude that the black universes un- 
der consideration are stable against axial perturbations 
in the static region. 

As we have shown, black universes at c = — 1 (6 = 
37rm/2) are stable against polar monopole perturbations. 
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FIG. 3: Effective potentiai as a function of tlie radial coordinate r (left figure) and of the tortoise coordinate r, (right figure) 
for axial gravitational perturbations of M-AdS wormholes with m = —1/3 for I = 2, c = 0.8 (no negative gap), 1.0, 1.8 (a 
negative gap), 2.2 (a deep negative gap). 




FIG. 5: Left panel: Effective potentials for radial perturbations of a black universe with I = Q, m = 1/3, c = —1 (red, top), 
c = —1.001 (green), c = —1.01 (blue, bottom). The potentials vanish at the horizon. Right panel: time-domain perturbation 
evolution for c = —1 (red, stable), and c = —1.001 (green, unstable). 
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Veff(r) 




In Table I, the asymptotic formula for the fundamental 
mode (the one that dominates at late times) is presented 
in units m = 1/3, so that b = 7r/2. 

The WKB formula (52) used here is developed for ef- 
fective potentials which have the form of a barrier with 
only one peak (see, e.g.. Fig. 9 for black universes), so 
that there are two turning points given by the equation 
— Veff = 0. Therefore, it cannot be used for effective 
potentials which have negative gaps, that is, for testing 
the stability of all questionable cases. 

VI. CONCLUSIONS 



FIG. 6: Effective potential for radial perturbations of a black 
universe with £ — 0, m = 1/3, c = —1.5 (red, top), c = —2 
(green), c — —3 (blue, bottom). The potentials vanish at the 
horizon. 

TABLE I: Quasinormal modes of axial perturbations of black 
universes with c — ~1 (m — 1/3). 



£ 


WKB6 


time-domain fit 


2 
3 
4 
5 
> 1 


1.168844 - 0.176095i 
1.590954 - 0.182539i 
1.997792 - 0.185591i 
0.191226(2£ + 1 - i) 


0.712299 - 0.157848i 
1.169144 - 0.176019i 
1.590956 - 0.182542i 
1.997801 - 0.185588i 



and their response to external perturbations is dominated 
at late times by the quasinormal (QN) frequencies. Sup- 
posing 4' oc e~'"*, quasinormal modes can be written in 
the form 

where a positive is proportional to the decay rate of 
a damped QN mode. The low-lying axial QN frequencies 
have the smallest decay rates in the spectrum and thus 
dominate in a signal at sufficiently late times. They can 
be calculated with the help of the WKB approach [38, 39]. 

Introducing Q = ui^ — Ves, the 6-th order WKB formula 
reads 

-^-Y.^'^-^+i n = 0,l,2..., (52) 

where the correction terms were obtained in [38, 39]. 
Here Qo and Qq are the value and the fc-th derivative of 
Q at its maximum with respect to the tortoise coordinate 
r*, and n labels the overtones. The WKB formula (52) 
was effectively used in a lot of papers (see, e.g., [40] and 
references therein). 

The WKB formula gives an accurate result for large 
multipole numbers (see Table I). Expanding the WKB 
formula it powers of £, we find the following asymptotical 
expression for c = — 1: 




We have developed a general formalism for analyzing 
axial gravitational perturbations of an arbitrary static, 
spherically symmetric solution to the Einstein-Maxwell- 
scalar equations where the scalar field, which can be both 
normal and phantom, is minimally coupled to gravity 
and possesses an arbitrary potential. This can be used 
for studying the stability and QNM modes of diverse 
charged and neutral scalar field configurations in general 
relativity. As to polar perturbations, we have restricted 
ourselves to the monopole mode, i.e., to spherically sym- 
metric (radial, for short) perturbations. 

We have applied this formalism to some electrically 
neutral wormholes and black holes supported by a phan- 
tom scalar field. The main results which were obtained 
here are: 

1. M-AdS wormholes described by the solution (3) 
with negative mass are shown to be unstable under 
radial perturbations, although the initial effective 
potential with a singularity at the throat is positive 
everywhere. These features are similar to those of 
anti-Fisher wormholes [14] which are twice asymp- 
totically flat. 

2. Black universes (i.e., regular black holes with de 
Sitter expansion far beyond the horizon) , described 
by the solution (3) without throats in the static re- 
gion, are shown to be stable only in the special 
case where the horizon coincides with the mini- 
mum of the area function R(r) (the parameters: 
c = — 1, b = 37rm/2) and unstable for c ^ —1. 

3. Quasinormal modes of axial perturbations have 
been calculated for stable black universes. 

Quite a lot of other problems of interest are yet to be 
studied. One can mention a full nonlinear analysis of per- 
turbations in all relevant cases. Next, the formalism de- 
veloped here for linear perturbations allows us to include 
the electromagnetic field into consideration and to study 
charged solutions with both normal and phantom scalar 
fields. Last but not least, using the well-known confor- 
mal mappings that relate Einstein and Jordan frames of 
scalar-tensor and curvature-nonlinear theories of gravity, 
one can extend the studies to solutions of these theories. 
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FIG. 7: Regularized effective potentials (left panel, top to bottom) and time-domain profiles (right panel, bottom to top) for 
radial perturbations of a black universe with £ = 0, m — 1/3, c — —0.1 (red), —0.2 (green), —0.3 (blue). 




FIG. 8; Regularized effective potentials (left panel) and time-domain profiles (right panel) for radial perturbations of a black 
universe with £ — 0, m — 1/3, c = —0.90 (blue), —0.95 (green), —0.99 (red). Larger negative values of c correspond to deeper 
negative gaps and later growing phase of the signal. 



fig 9 




FIG. 9: The effective potential (left figure) and the time-domain profile (right figure) for axial gravitational perturbations of 
black universes for £ — 2, m = 1/3, c = —1. The potential vanishes at the horizon r = 0. 
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